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(57) Abstract: An image of an object is synergistically reconstructed using two or multiple imaging modalities. A first reconstructed 
image, showing structural information of the object is produced using a first imaging modality. The first reconstructed image is 
segmented, and known optical properties of the object are then mapped to the first reconstructed image. Optical signal emissions 
from the object are detected and registered with the first reconstructed image. A second reconstructed image volume is then produced 
using a second imaging modality, based on the mapped optical properties after registration between the first image and the data from 
the second modality. The second reconstructed image depicts some optical property, such as a bioluminescent source distribution, 
or optical properties, such as, attenuation and scattering properties, of the object. 
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BACKGROUND 

This invention relates to multi-modality-based systems and methods for 
detecting an optical property distribution, such as a light-emitting source distribution, in 
1 0 multiple dimensions as well as systems and methods for reconstructing such an image 
from the detected signals from the distribution based on data from a tomographic 
imaging modality, including but not limited to computed tomography (CT) or micro- 
CT. 

There are many "emission-detection 95 imaging techniques known in the art, such 
15 as bioluminescent imaging. However, such current imaging techniques are limited to 
the projective imaging mode or external excitation of the internal light source through 
external energy sources along selected paths. Therefore, three-dimensional structures 
and localization of an internally derived light source, such as one not reliant on external 
energy excitation, cannot be resolved with high quantitative accuracy both in terms of 
20 spatial location and localized activity. 

Diffuse computed tomography (CT), another known imaging technique, 
computes distributions of absorption and scattering coefficients from scattered light 
through an object. Typically, intensity-modulated light sources are used. It is well 
known that diffuse CT will generally produce low image resolution, particularly as 
25 background heterogeneity increases. 

It would therefore be desirable to combine an optical imaging technique, such 
as a light emission technique, specifically bioluminescent imaging, with a scanning 
technique that allows the evaluation of two and three dimensional structural 
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information, such as computed tomography scanning or magnetic resonance imaging, 
to produce a reconstructed image having better image resolution. 

SUMMARY 

The present invention is directed to multi-modal imaging systems and methods 
5 that reconstruct images via fundamental and synergistic utilization of multi-model data. 
According to an exemplary embodiment, an image volume may be reconstructed in a 
first tomographic modality, optical properties from a database may be mapped to the 
image volume, and the image may then be reconstructed tomographically in another 
modality, based on the optical properties. 
1 0 According to one embodiment, bioluminescent CT and CT/micro-CT 

combinations may be used, but other system configurations are possible. Some 
embodiments may include a magnetic resonance imaging (MRI) scanner or micro -MRI 
scanner in conjunction with a fluorescent tomographic scanner. The imaging 
techniques and algorithms described herein are exemplary only, and other methods of 
1 5 combining data from two or more tomographic scanners may be used. 

Some embodiments may be capable of various resolutions depending on 
scanning times, possess extremely high photon detection sensitivity for mapping gene 
expression, and/or embody hardware and/or software technology for image 
reconstruction, registration, and analysis. Some embodiments may have the advantage 
20 of being configured to rapidly collect data with a higher signal-to-noise ratio. 

According to one embodiment, bioluminescent imaging may be rendered in a 
two-or three-dimensional tomographic modality. In embodiments directed to 
bioluminescence, emitted photons can be collected from multiple three-dimensional 
directions with respect to an animal marked by bioluminescent compounds including 
25 reporter luciferases. 

According to some embodiments, a CT or micro-CT scanner may be integrated 
with a bioluminescent imaging system. The bioluminescent imaging system may also 
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be combined with other imaging systems which provide information regarding the 
distribution of tissue structures in vivo, in situ, or ex vivo. 

In alternative embodiments, an object may be serially scanned using each 
modality in turn. In still further embodiments, the object may be transported between 
5 scanning modalities. Optionally, one or more registration marks may be placed on the 
object to coordinate positions between scanning modalities. The surface of the object 
may also be optically reconstructed for the registration purpose. 

In some embodiments, information associated with x-ray CT imaging and 
bioluminescent imaging may be used together to estimate light scatter and/or other 

1 0 optical properties of tissue and thereby reconstruct a three-dimensional emission image 
volume registered to corresponding CT or micro-CT imaging of anatomical and 
pathological structures. As non-limiting examples, the system may be used to generate 
images of structures, such as bioluminescent sources, lungs and various tumors. 

According to some embodiments, intra-organ localization of gene transcription 

15 activity may be performed with resolution capable of differentiating, for instance, gene 
expression in the central pulmonary airways (out to approximately the 5th-7th 
generation) versus parenchymal activity. Also, localization of parenchymal activity in 
terms of sub-lobar regions may be performed. As a non-limiting example, small animal 
imaging, in particular mouse imaging, may be performed. In other examples, the 

20 systems and methods may be used for other biomedical applications where 

bioluminescent signals are detectable. Some embodiments are especially suited for 
small animal imaging at molecular levels. For example, genetic activity in a particular 
organ system may be imaged. 

By integrating x-ray and optical imaging, better optical tomography image 

25 quality can be achieved that would not be possible with a stand-alone optical system. 
From a corresponding x-ray CT image volume or image volume generated by other 
imaging energy sources, knowledge of the underlying distribution of optical scatters 
can be derived. This information is useful in reconstruction of images from optical 
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data. Specifically, emitting source distributions may be directly solved for, obviating 
the need for reconstruction of optical properties in three dimensions. 

According to exemplary embodiments, the combined use of x-ray CT and 
BLCT transforms the nonlinear optical CT problem into an easier linear problem. 
5 Therefore, the reconstruction of image data from a bioluminescent CT scanner may be 
significantly improved. 

One embodiment includes a system processor that supports the desired 
functionality as described in detail below and a system data store (SDS) that stores data 
associated with the needed functionalities, such as image data and reconstruction. The 
10 system processor may be in communication with the SDS via any suitable 
communication channel(s). 

The SDS may include multiple physical and/or logical data stores for storing the 
various types of information used. Data storage and retrieval functionality can be 
provided by either the system processor or one or more data storage processors 
15 associated with the SDS. The system processor may include one or more processing 
elements that are adapted or programmed to support the desired image storage, 
reconstruction and/or other functionality. 

Accordingly, one method of image reconstruction includes a variety of steps 
that may, in certain embodiments, be executed by the environment summarized above 
20 and more fully described below or be stored as computer executable instructions in 

and/or on any suitable combination of computer-readable media. The steps can include 
but are not limited to performing tomographic reconstruction of an image volume in 
one modality, mapping optical properties to that volume from a database, and 
performing tomographic reconstruction in another modality based on the mapped 
25 optical properties. 

Additional advantages of the invention will be set forth in part in the description 
which follows, and in part will be obvious from the description, or may be learned by 
practice of the invention. It is to be understood that both the foregoing general 
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description and the following detailed description are exemplary and explanatory only 
and are not restrictive of the invention as claimed. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 The accompanying drawings, which are incorporated in and constitute a part of 

this specification, illustrate embodiments of the invention and together with the 
description, serve to explain the principles of the invention. 

FIG. 1 depicts an exemplary bioluminescent imaging device (left panel) with an 
anatomic imaging device (micro CT scanner shown in the right panel). 
1 0 FIG. 2 depicts an exemplary single source bioluminescent CT scanner. The 

rotating stage and light sensitive camera is shown on the bottom shelf of the cart. The 
side of the light tight enclosure has been removed for system visualization. 

FIG. 3 depicts bioluminescence emitted from the lungs of a mouse following 
exposure to an adenovirus which has delivered the primary gene and reporter gene 
15 (producing luciferase) to the lungs. The emitted light can be seen from multiple angles 
according to exemplary embodiments. This is a precursor to being able to reconstruct 
the 3D distribution of the light source (s). 

FIG. 4 depicts an exemplary micro-CT image of a lung showing structural 
components of the mouse thorax with resolution down to the alveolar level. This sort 
20 of anatomic image serves to provide the knowledge of the distribution of light 
scatterers. 

DETAILED DESCRIPTION 

One or more exemplary embodiments are now described in detail hereinbelow 
25 and in the attachments hereto. Referring to the drawings, like numbers indicate like 
parts throughout the views. As used in the description herein and attachments hereto, 
the meaning of "a," "an," and "the" includes plural reference unless the context clearly 
dictates otherwise. Also, as used in the description herein and attachments hereto, the 
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meaning of "in" includes "in" and "on" unless the context clearly dictates otherwise. 
Finally, as used in the description herein and attachments hereto, the meanings of "and" 
and "or" include both the conjunctive and disjunctive and may be used interchangeably 
unless the context clearly dictates otherwise. 
5 Ranges may be expressed herein as from "about" one particular value, and/or to 

"about" another particular value. When such a range is expressed, another embodiment 
includes from the one particular value and/or to the other particular value. Similarly, 
when values are expressed as approximations, by use of the antecedent "about," it will 
be understood that the particular value forms another embodiment. It will be further 
understood that the endpoints of each of the ranges are significant both in relation to the 
other endpoint, and independently of the other endpoint. 

The present invention relates to systems and methods for detecting a light- 
emitting source distribution as well as systems and methods for reconstructing an image 
from the detected light signals and tomographic images obtained from one or more 
other modalities, such as an image volume from CT or micro-CT. Some embodiments 
may include one or more cameras arranged, e.g., symmetrically, on a spherical surface 
to detect a light emitting source distribution in three dimensions. Alternative 
embodiments may include asymmetrical camera arrangements and/or other three- 
dimensional surface arrangements. In some embodiments, other optical mechanisms 
can be used to intercept and direct signals to the cameras including, but not limited to, 
mirror and/or fiber systems. 

Some further embodiments may detect and record bioluminescent emissions 
and/or fluorescent emissions. This image data, along with associated x-ray CT images 
of the same object, can be used to reconstruct a three-dimensional emission image 
volume and register the bioluminescent CT image to a corresponding x-ray CT or 
micro-CT image volume of anatomical and pathological structures. In some such 
embodiments, the bioluminescent (or fluorescent) reconstruction process can be 
enhanced through the use of knowledge gained from x-ray CT or other anatomic 
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infoimatibn^ of other imaging devices including, but not limited to, 

MRI or ultrasound. As a non-limiting example, emitted photons can be collected from 
multiple directions in three dimensions with respect to a living animal or any other light 
emitting structure of interest marked by bioluminescent reporter luciferases or 
5 fluorescent sources. In some embodiments, a lung and/or various tumors may be 
imaged. 

According to exemplary embodiments, optical properties of an object are 
presumed to be already known, and the photon-emitting source distribution may be 
computed based on these known optical properties. Therefore, the imaging models for 

1 0 the systems and methods according to exemplary embodiments are approximately 

linear, while that for the conventional techniques, such as diffuse CT, are nonlinear and 
correspondingly more difficult to solve. 

By combining a system for detecting light emission from multiple angles of 
view simultaneously or sequentially with an imaging modality which allows the 

1 5 evaluation of two and three dimensional structural information, such as micro x-ray CT, 
the anatomic and/or structural details gained from the micro x-ray CT can be used to 
estimate the distribution of light scattering structures for purposes of directing the 
computed tomographic calculations required to create BLCT cross-sectional or 
volumetric images. Such a system may enable, as non-limiting examples, both the 

20 calculation of the computed tomograms of chemo- luminescence and the linking of the 
computed tomograms of chemo-luminescence to the highly detailed anatomic image 
sets derived from the x-ray CT imaging, hi some embodiments, the tomographic 
reconstruction of bioluminescence can provide important added detail regarding 
regional location of reporter gene activity. By knowing the location of reporter gene 

25 activity and having micro-resolution images of anatomy, a user can follow the link 
between gene activation and pathologic processes. 
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Typical Storage and Processing Architecture 

In one exemplary embodiment, the imaging and reconstruction system includes 
a system processor potentially including multiple processing elements. The term 
processing element may refer to (1) a process running on a particular piece, or across 
5 particular pieces, of processing hardware, (2) a particular piece of processing hardware, 
or either (1) or (2) as the context allows. Each processing element can be supported via 
a standard general purpose processor such as an Intel-compatible processor platforms 
preferably using at least one CELERON, PENTIUM, XEON, ITANIUM (Intel Coip., 
Santa Clara, CA) class processor; alternative processors such as MIPS (MIPS 

10 Technologies, Mountain View, CA) or UltraSPARC (Sun Microsystems, Palo Alto, 
CA) could be used in other embodiments. The system processor, or the one or more 
processing elements thereof, can include one or more field programmable gate arrays 
(FPGAs), programmable digital signal processors (DSPs) and/or application specific 
integrated circuits (ASICs) configured to perform at least a portion of the functionality 

1 5 according to the present invention. In other embodiments, an embedded 

microprocessor can be used such as, but not limited to 5 an ARM (ARM, Carlsbad, CA) 
processor core. 

In some embodiments, the system processor can include a combination of 
general purpose processors, ASICs, DSPs and/or FPGAs. In some embodiments, the 
20 systems and methods of the present invention, as described above, can be distributed 
across multiple processing elements. In some such embodiments, aspects of the 
functionality or portions thereof may be executed in series or in parallel; particular 
functionality or portions thereof executed a multiplicity of times may also occur in 
series or parallel. 

25 In a system processor including at least one general purpose processor, the 

general purpose processor typically runs an appropriate operating system such as 
WINDOWS/NT, WINDOWS 2000 or WINDOWS/XP (Microsoft, Redmond, WA), 
IRIX (Silicon Graphics, Mountain View, CA), SOLARIS (Sun Microsystems, Palo 
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Alto, CA), or LINUX (or other UNIX variant). In one embodiment, the Windows 2000 
operating system is used. 

The SDS may include a variety of primary and secondary storage elements. In 
one embodiment, the SDS can include random access memory (RAM) as part of the 
5 primary storage; the amount of RAM might range from 512 MB to 4 GB in some 

embodiments. The primary storage can, in some embodiments, include other forms of 
memory such as cache memory, registers, non-volatile memory (e.g., FLASH, ROM, 
EPROM, etc.), etc. 

The SDS can also include secondary storage including single, multiple and/or 

10 varied servers and storage elements. For example, the SDS can use internal storage 

devices connected to the system processor. In embodiments where a single processing 
element supports all of the system functionality, a local hard disk drive can serve as the 
secondary storage of the SDS, and a disk operating system executing on such a single 
processing element can act as a data server receiving and servicing data requests. A 

15 system bus can serve as the communication channel between the system processor and 
the SDS (typically, at least RAM and the hard disk drive). 

It will be understood by those skilled in the art that the different information 
used in the imaging and image reconstruction processes and systems according to the 
present invention can be logically or physically segregated within a single device 

20 serving as secondary storage for the SDS; multiple related data stores accessible 

through a unified management system, which together serve as the SDS; or multiple 
independent data stores individually accessible through disparate management systems, 
which may in some embodiments be collectively viewed as the SDS. The various 
storage elements that comprise the physical architecture of the SDS may be centrally 

25 located or distributed across a variety of diverse locations. 

The architecture of the secondary storage of the system data store may vary 
significantly in different embodiments. In several embodiments, database(s) are used 
to store and manipulate the data; in some such embodiments, one or more relational 
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database management systems, such as DB2 (IBM, White Plains, NY), SQL Server 
(Microsoft, Redmond, WA), ACCESS (Microsoft, Redmond, WA), ORACLE (Oracle 
Corp., Redwood Shores, CA), Ingres (Computer Associates, Islandia, NY), MySQL 
(MySQL AB, Sweden) or Adaptive Server Enterprise (Sybase Inc., Emeryville, CA), 
5 may be used in connection with a variety of storage devices/file servers that may 

include one or more standard magnetic and/or optical disk drives using any appropriate 
interface including, without limitation, ATA, IDE and SCSI. In some embodiments, a 
tape library such as available from Exabyte Corporation (Boulder, CO), a storage 
attached network (SAN) solution such as available from EMC, Inc. (Hopkinton, MA), a 

1 0 network attached storage (NAS) solution such as available from Network Appliances 
(Sunnyvale, CA), or combinations thereof may be used. In other embodiments, the 
data store may use database systems with other architectures such as object-oriented, 
spatial, object-relational or hierarchical. 

Instead of, or in addition to, those organization approaches discussed above, 

1 5 certain embodiments may use other storage implementations such as hash tables or flat 
files or combinations of such architectures. Such alternative approaches may use data 
servers other than database management systems such as a hash table look-up server, 
procedure and/or process and/or a flat file retrieval server, procedure and/or process. 
Further, the SDS may use a combination of any of such approaches in organizing its 

20 secondary storage architecture. 

The SDS communicates with the system processor by one or more 
communication channels. Multiple channels can be involved in some embodiments for 
supporting communication between processing elements of the system processor and 
portions of the SDS. Such channels can include without limitation computer network, 

25 direct dial-up connection, dedicated connection, direct or indirect connection such as 
via a bus connection, parallel or serial connection, USB connection, null modem 
connection or wireless connection utilizing an appropriate communication protocol 
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such as BLUETOOTH, IRDA, 802. 1 lb or other suitable channel as would be known to 
those skilled in the art. 

All forms of data, including raw, intermediate, and computed can be stored on 
one or more SDS either temporarily or permanently. In particular, the SDS can store, 
5 without limitation, image data, including volumetric image data, reconstruction 
intermediate data, final reconstructed imaging data, imaging parameters, and 
reconstruction parameters. Further, the SDS may, in some embodiments, store 
instructions for performing the various imaging and reconstruction tasks, or portions of 
such tasks. 

10 

Light Sensitive Cameras 

In one embodiment, ten CCD cameras can be arranged at the center of each 
identical face of a dodecahedron, except for the two facing the front and back ends of 
an object to be imaged. Each of the 10 cameras can point to the iso-center where the 

15 object can be fixed on a holder. The imaging geometry can be implemented by a 
structure holding each camera in its fixed position. Data from the cameras can be 
transmitted to one or more processing elements for further processing and image 
reconstruction. A light-free housing can be used to house the imaging cameras. Once 
the frame of the imaging device is arranged, a camera can be mounted at any spot of the 

20 12 nominal positions on the dodecahedron. One or more cameras can be geometrically 
and photographically calibrated with reference phantoms. An optical surface scanner 
can be combined with the imaging frame for the registration purpose. 

One skilled in the art will recognize that other arrangements of cameras are 
possible, including geometrically symmetrical and asymmetrical configurations, with or 

25 without appropriate optical paths such as obtained through use of mirrors or fiberoptic 
relay paths. In some embodiments, one or more cameras can be rotated around an 
object of interest. Alternatively, or in combination, the subject of the imaging can be 
rotated on one or more axis. As a non-limiting example, such embodiments can be 
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used in cases where light emission is unstable. In this case, dynamic bioluminescent 
and/or fluorescent tomographic imaging would be feasible. 

In one embodiment, one or more cooled back-thinned integrating CCD cameras 
can be used for imaging. The camera package can include a 2.2L or other appropriately 
5 sized end-on liquid nitrogen dewar for cooling. Alternatively, an omni-directional 

dewar can be used to allow cameras to be mounted in any orientation while keeping the 
dewar right-side up. Li some embodiments, the imagers can be sensitive to one or more 
bioluminescent sources of different spectral characteristics. Other types of light 
sensitive cameras can also be used. Analog films can be used with appropriate 

10 manipulations. 

A living organism, or other structure of interest, can be scanned using a multi- 
detector spiral CT, another appropriate method of imaging known to one skilled in the 
art, or a micro x-ray CT scanner. From this imaging, a distribution of optical properties 
of the object can be derived to guide associated bioluminescent and/or fluorescent 

1 5 tomographic imaging. 

Camera Control 

According to exemplary embodiments, one or more camera control elements 
may be used. A camera control element can include one or more processing elements 

20 and can be in communication with the SDS. The imaging cameras of the present 
invention can be in communication with the one or more camera control elements. 
Camera control elements can perform digitization of output from cameras and other 
processing as appropriate. Relevant imaging parameters can be controlled by the 
camera control elements. As non-limiting examples, imaging parameters such as focus, 

25 exposure time, aperture, can be configured. Additional parameters known to one 
skilled in the art can also be configured as appropriate. 

In some embodiments, all cameras can be controlled by a single controller 
element. In other embodiments, a single controller element may be used to control 



WO 2004/081865 



PCT/US2004/006360 



13 

multiple cameras. In such an embodiment, cameras may be arbitrarily grouped into 
arrays and each array can be controlled by a single controller element. Each camera 
acquisition chain can then operate independently. Still further embodiments can 
include a master controller element to provide control and synchronization for an entire 
5 camera array. In one embodiment, one camera is controlled by a single controller 

element. One skilled in the art will recognize that other configurations of cameras and 
controller elements are possible. 

Some embodiments can include one or more hierarchies of control elements 
and/or processing elements. Different levels of hierarchy can perform the same or 
10 different functions. In one embodiment, images from low-level camera control 

elements are passed to higher level processing elements and/or control elements for 
additional processing and/or image reconstruction, as further described below. Control 
elements and processing elements can be communicatively coupled via any suitable 
communication means including computer network, wired or wireless direct link, bus 
1 5 connection and as further described below. 

In one embodiment, individual cameras can be controlled with an external 
camera control element in communication with a PCI controller card in communication 
with one or more processing elements. The PCI controller card and the camera control 
element can be communicatively coupled via an RS-422 link or other suitable serial or 
parallel connection. The external camera control element and the camera can be 
communicatively coupled via a parallel interface. In one such embodiment, cameras 
output analog signals which are digitized by one or more external camera control 
elements. To minimize noise, the length of the link between the camera and the 
controller box can be minimized. 

In one embodiment, the light-emitting source data acquisition process can 
include one or more of the following steps: (1) reset and/or initialize cameras; (2) 
execute a programmable, configurable, or manual shutter open on one or more cameras; 
(3) execute a programmable, configurable, or manual shutter close on one or more 
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cameras; (4) transfer image data to one or more camera control elements, processing 
element, or SDS; and (5) store the images in the SDS. Additional embodiments can 
have shutter times configured to occur simultaneously across multiple cameras and/or 
automatically after a predetermined time delay. 

5 

Image Reconstruction Using an Iterative Method or Other Methods 

Given the ill-posed nature of the imaging and sampling geometry, an iterative 
image reconstruction approach may be used utilizing prior knowledge on the 
distribution to be reconstructed. The iterative approach can be used in the case of 
10 incomplete and/or noisy data. Also, the iterative approach easily accommodates prior 
knowledge and imaging physics. 

An interface from one or more camera control elements and/or processing 
elements to one or more reconstruction engines can be provided. The systems of the 
present invention can acquire optical properties of the object being imaged and then 
15 compute the light-emitting source distribution based on the optical properties. 

Therefore, the imaging model is approximately linear. In one embodiment, x-ray CT 
data is used to regularize the BLCT reconstruction problem and transform it from a 
nonlinear one to linear one and thereby greatly simplify it. 

Even with attenuation and scattering taken into account based on a CT or 
20 micro-CT image volume, a discrete BLCT imaging model can still be linearly 

expressed as Ax = b , where the observed data b = (b l , • • • , b M ) e R M , original emitting 
source distribution x = (x 1 , • ■ • , x M ) s R N , and a known non-zero M x N matrix 
A = {Aj ) • The coefficients of the matrix A depend on the anatomical structures and 

their optical properties according to the classic Radiative Transfer Equation (using the 
25 Monte Carlo method) or the well-known diffusion approximation or another 

appropriate method. The systems and methods according to exemplary embodiments 
can reconstruct the image x from the data b . 
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A generalized BLCT algorithm according to one embodiment cair include one 
or more of the following steps: (1) reconstruction and segmentation of an x-ray CT 
image volume, (2) association of optical properties to each segmented region in the x- 
ray CT volume based on a library of optical properties, (3) determination of coefficients 
5 of the forward imaging matrix A = {A iy ) based, e.g., on Monte Carlo simulations, (4) 

reconstruction of the emitting source distribution x by inverting the matrix A, subject 
to the constraints imposed by the segmented anatomical structures, their properties, and 
known features of the underlying source distribution (such as the homogeneity or 
parametric form of the source distribution, the shape or intensity of the source). 

1 0 According to exemplary embodiments, bioluminescent emissions, and well as other 
light-emitting source distributions, such as fluorescent source distributions, may be 
detected and reconstructed in this way or another alternative way. 

As non-limiting examples, the optical properties of step two can include 
absorption coefficients, scattering coefficients, scattering anisotropy, indices of 

15 refraction, and other appropriate parameters known to one skilled in the art. Monte 
Carlo simulations can be used to predict bioluminescent signals and construct the 
matrix A based on the CT/micro-CT image volume of the object. Other methods, such 
as finite element methods and meshfree methods, may be also used for this purpose. 
After image segmentation, optical properties can be assigned to each segment based on 

20 a library of optical properties. 

In one embodiment, both the ordered-subset expectation maximization (OS- 
EM) and the ordered-subset version of the simultaneous algebraic reconstruction 
technique (OS-SART) schemes for BLCT can be implemented. A roughness penalty 
method for BLCT or other method known to one skilled in the art can also be used. 

25 Although an iterative method may be most suitable to the image reconstruction 

task in one embodiment, other image reconstruction methods can be used. Even 
further, the iterative procedure described above is only an example, and should not be 
interpreted as a limiting description. 
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As far as image reconstruction methods are concerned, it is emphasized that 
there are multiple options or possibilities. In addition to an iterative reconstruction 
strategy as described above, numerical solutions to the Radiative Transfer Equation or 
its approximation, such as the diffusion equation, can be useful as well. A fast analytic 
5 method would be very useful in practice. In one embodiment, an analytic approach 
known as the Kirchhoff approximation may be adapted for bioluminescent tomography 
of diffuse media with an embedded source distribution. Other numerical methods, such 
as finite element methods and/or meshfree methods, are also feasible for the same 
purpose. 

10 

CT/Micro-CT Scanner 

Any state-of-the-art micro-CT scanner can be used according to exemplary 
embodiments. In some embodiments, the hnTek Micro CAT II described, e.g., at 
http://www.imtekinc.coni/litml/niicrocat ii specifications.htmL or the SkyScan-1076 

15 in-vivo micro-CT system described, e.g., at 

http://www.skvscan.be/next/spec 1076.htm . can be used. 

In other embodiments 2 the data acquisition system can include one or more of 
the following items: a dedicated embedded data acquisition and control computer, two 
130 kVp ultra-high resolution jnfocus X-ray systems, two 100/50 mm dual-field image 

20 intensifiers, and two 2048 x 2048 CCD cameras. The scanner can include one or more 
processing elements in communication with the SDS, a multi-axis precision scanner 
and specimen manipulator with linear servo drives, remotely configurable motorized 
source-detector geometries, a signal and power slip ring for continuous rotation, and a 
means to move data between acquisition and processing. The slip ring can have two 

25 independent capacitively coupled data transmission channels with full-duplex fibre 
channel interfaces. The one or more processing elements, the slip ring data channels 
and the SDS can be communicatively coupled. In one embodiment, they can be 
connected via one or more fiber optic cables. 
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One embodiment may be built on an optical grade table for vibration isolation 
and precision alignment. An imaging chain of source arrays, detector arrays, and 
accessories, alone or in combination, can be mounted on a rotating plate which is in 
turn supported by an open bearing and rigid stand. The axis of rotation can set to any 
5 appropriate angle including vertical and horizontal. The geometry of each imaging 

chain can be individually configured to suit a wide variety of operating modes. Each x- 
ray tube and image intensifier can be moved radially, while each image intensifier can 
also be moved laterally. The object can rest horizontally in a holder mounted to a 
linear axis with a certain amount of axial travel room for slice positioning. One 
10 exemplary embodiment is capable of achieving spatial resolution of 100 lp/mm for 
excised samples, and temporal resolution of 1.8 seconds for objects up to 120 mm in 
diameter. The system can be configured to allow a wide range of intermediate 
combinations of scan time and spatial resolution. 

1 5 System Integration 

The light-emitting source distribution CT device and the anatomic imaging 
scanner, such as a micro x-ray CT scanner, can be electronically and mechanically 
integrated but need not be in all embodiments. In one embodiment, the hardware 
structures of the two imaging units can share a table and/or a holder attached to a table. 

20 This embodiment can allow the translation of an object for x-ray CT scanning to be 

extended into the light-emitting CT device in a precise and/or repeatable fashion. Some 
embodiments may be configured to optimize and integrate software packages for 
Monte Carlo simulation (another kind of simulation, such as that based on finite 
element computation), CT and/or micro-CT data preprocessing and reconstruction, 

25 BLCT reconstruction, image visualization and analysis. A user interface to perform 

and/or to configure such functions can also be provided in some embodiments; in some 
such embodiments, the user interface can further allow viewing of results and may 
allow control of parameters with respect to such viewing. Any software capable of 
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performing such functions can be implemented on one or more processing elements. 

Exemplary Applications 

The following applications are intended as illustrative examples only and are 
5 not limiting of the invention. According to exemplary embodiments, advanced 
imaging, such as lung imaging, is enabled in that the structural and function 
information can be obtained concurrently at the molecular level, and can be evaluated 
on a regional, sub-lobar basis. This combination allows simultaneous examination of 
gene expression and anatomic structures and improves understanding of the human 
10 lungs. 

Exemplary embodiments may be used in gene therapy imaging, to probe the 
distribution of the administered gene, reporter genes, such as those producing 
luciferase, can be included in the transfecting virus. These genes cause the emission of 
light, enabling the functional gene to be identified within the target tissue. 

15 Exemplary embodiments may also be useful in evaluating transgene expression 

in the lung; gene transfer vectors; gene transfer to the respiratory epithelium of mice; in 
vivo bioluminescence imaging; understanding the site of transgene expression in the 
lung; human lung lobe imaging and sheep-based emphysema model evaluation; 
understanding the site of gene therapy, and its consequences; and understanding the 

20 pathophysiology of airway vs. alveolar infection. 

The embodiments described above are given as illustrative examples only. It 
will be readily appreciated by those skilled in the art that many deviations and other 
applications may be made from the specific embodiments disclosed in this specification 
without departing from the scope of the invention. 
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What is claimed is: 

1 . A method for reconstructing an image of an object, comprising: 
imaging the object using a first imaging modality to produce a first 

reconstructed image; 

mapping optical properties of the object to the reconstructed image volume; and 
detecting optical signals emitted from the object using a second imaging 

modality to produce a second reconstructed image, based on the mapped optical 

properties. 

2. The method of claim 1 5 wherein the first reconstructed image shows two or 
three dimensional structural details of the object. 

3. The method of claim 1, wherein the second reconstructed image shows two or 
three dimensional distribution of light emission from the object. 

4. The method of claim 1 5 wherein the second reconstructed image is reconstructed 
for multiple types of source distributions with various spectral characteristics. 

5. The method of claim 1, wherein the second reconstructed image is reconstructed 
from a single or multiple angles of view. 

6. The method of claim 1, wherein the second reconstructed image is reconstructed 
using an iterative or analytic approach. 

7. The method of claim 1, wherein the step of detecting optical signals uses 
sensors. 
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8. The method of claim 7, wherein the step of detecting optical signals also uses 
optical path components. 

9. The method of claim 1, wherein the second reconstructed image shows cross- 
sectional or volumetric views of the object or quantitative features of underlying source 
distributions of the object. 

10. The method of claim 1, wherein the optical properties include at least one of 
absorption coefficients, scattering coefficients, scattering anisotropy, indices of 
refraction, and features of underlying sources. 

1 1 . The method of claim 1 5 wherein the first imaging modality includes at least one 
of x-ray computed tomography, micro computed tomography, magnetic resonance 
imaging, and ultrasound. 

12. The method of claim 1 9 wherein the second imaging modality includes at least 
one of bioluminescent tomography and fluorescent tomography. 

13. The method of claim 1, further comprising segmenting the first reconstructed 
image into regions, wherein the step of mapping maps the optical properties to each 
segmented region of the image. 

14. The method of claim 1, further comprising registering the first reconstructed 
image with the detected optical signals before producing the second reconstructed 
image. 

15. The method of claim 14, wherein the step of registering uses a landmark-based 
method, a land-mark free method, or an optical surface imager method. 
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16. A system for reconstructing an image of an object, comprising: 

a first imaging device for imaging the object using a first imaging modality to 
produce a first reconstructed image; 

a library of optical properties of the object; 

a processor for mapping the optical properties of the object to the first 
reconstructed image; and 

a second imaging device for detecting optical signals emitted from the object 
using a second imaging modality to produce a second reconstructed image, based on 
the mapped optical properties. 

17. The system of claim 16, wherein the first reconstructed image shows two or 
three dimensional structural details of the object. 

18. The system of claim 16, wherein the second reconstructed image shows two or 
three dimensional distribution of light emission from the object. 

19. The system of claim 16, wherein the second reconstructed image is 
reconstructed for multiple types of source distributions with various spectral 
characteristics. 

20. The system of claim 16, wherein the second reconstructed image is 
reconstructed from a single or multiple angles of view. 

21 . The system of claim 16, wherein the second reconstructed image is 
reconstructed using an iterative or analytic approach. 
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22. The system of claim 16, wherein the second imaging device uses sensors for 
detecting the optical signal emissions. 

23. The system of claim 22, wherein the second imaging device further comprises 
optical path components. 

24. The system of claim 16, wherein the second reconstructed image shows cross- 
sectional or volumetric views or quantitative features of the underlying source 
distribution(s). 



25. The system of claim 16, wherein the optical properties include at least one of 
absorption coefficients, scattering coefficients, scattering anisotropy, indices of 
refraction, and features of underlying sources. 

26. The system of claim 16, wherein the first imaging modality includes at least one 

of x-ray computed tomography scan, micro computed tomography scan, magnetic 
resonance imaging, and ultrasound. 



27. The system of claim 16, wherein the second imaging modality includes at least 
one of bioluminescent tomography and fluorescent tomography. 

28. The system of claim 16, wherein the processor segments the first reconstructed 
image into regions and maps the optical properties to each segmented region of the 
image. 



29. The system of claim 16, wherein the processor registers the first reconstructed 
image with the detected optical signals before the second reconstructed image is 
produced. 
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30. The system of claim 29, wherein the processor performs registration using a 
landmark-based method, a landmark-free method, or an optical surface imager based 
method. 
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Figure 2. 
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Figure 3 
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Figure 4 



